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ABSTRACT 

We analyze the quasar two-point correlation function (2pCF) within the redshift in- 
terval 0.8 < z < 2.2 using a sample of 52303 quasars selected from the recent 7th 
Data Release of the Sloan Digital Sky Survey. Our approach to 2pCF uses a concept 
of locally Lorentz (Fermi) frame for determination of the distance between objects and 
permutation method of the random catalogue generation. Assuming the spatially flat 
cosmological model with given rt\ — 0.726, we found that the real-space 2pCF is fitted 
well with the power- low model within the distance range 1 < a < 35 Mpc with the 
correlation length tq = 5.85±0.33 h'^ Mpc and the slope 7 = 1.87±0.07. The redshift- 
space 2pCF is approximated with sq = 6.43 ± 0.63 Mpc and 7 = 1.21 ± 0.24 for 
1 < s < 10 /i-i Mpc, and sq = 7.37 ± 0.81 h'^ Mpc and 7 = 1.90 ± 0.24 for 
10 < s < 35 Mpc. For distances s > 10 Mpc the parameter describing the 
large-scale infall to density inhomogeneities is /3 — 0.63 ± 0.10 with the linear bias 
b = 1.44 ± 0.22 that marginally (within 2a) agrees with the linear theory of cosmo- 
logical perturbations. We discuss possibilities to obtain a statistical estimate of the 
random component of quasars velocities (different from the large-scale infall) . We note 
rather slight dependence of quasars velocity dispersion upon the 2pCF parameters in 
the region r < 2 Mpc. 
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1 INTRODUCTION 

Spatial distribution of quasars is one of the main sources 
of observational information about the largescale structure 
of the Universe; this is at present almost the only source 
of statistical information about matter inhomogeneity and 
clustering at cosmological redshifts. The most important re- 
sults on quasars clustering were obtained using two largest 
quasar surveys up to date: the 2-degree Field Quasar Survey 
with 2QZ cat a logue as a result (http://www.2dfquasar.org, 
ICroom et all (| 19981 )') and the Sloan Digital Sky survey 
(SDSS; http://www.sd ss.org), which has been completed 
with the 7th release jAbazadiian et all 20091'). The two- 
point correlation functions (2pCF) (Peebles" 1980') of quasars 
5 (r) are important characteristics of matter spatial inhomo- 
geneity that may be compared to cosmological theories of 
structure formation. Development of 2dF and SDSS surveys 
gave a powerful incentive to investigation of v arious aspects 
of correlation functions of quasars (see, e.g., ICroom et al.l 
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(|2005l ): lMvers" et all (|2006l ) and references therein). In the 
present paper we use the 7-th data release of SDSS to study 
2pCF of quasars. 

The main problems concerning reconstruction of the 
cosmological mass distribution from redshift surveys are as 
follows. The first problem is that the surveys of extragalactic 
objects give us an information only about distribution of the 
luminous matter wh ich is biased relative to the dark matter 
(|Dekel fc Reed 1 19871 ). Biasing may depend on the physical 
peculiarities of extragalac ti c obj e cts, e.g., on morpholog- 
ical type (|Einasto et all 12007 : IRoss. Brunner fc Mverj 
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of linear biasing means that the density variations of certain 
kind of objects (in our case quasars) is proportional to that 
of the whole matter; therefore 



C(r) = &^Cm(0, 



(1) 
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where ^m{r) is 2pCF of the whole matter, ^(r) is 2pCF of 
quasars and b is the bias parameter. 

The second problem of the 3-dimensional analysis of 
the matter distribution is due to the fact that the observed 
redshifts of extragalactic objects are 'contaminated' by mea- 
surement errors and non-Hubble motions. The distances cal- 
culated from these redshifts without taking into account the 
unknown peculiar velocities are called distances in redshift- 
space (in contrast to the real space). As our Universe is 
isotropic, the correlation function must be spherically sym- 
metric in the real-space. But in the redshift-space it appears 
to be distorted. On smaller scales the profile of galaxies 
2pCF is stretched along the line of sight ('Finger of God' ef- 
fect) due to virial velocities of objects inside the galaxy clus- 
ters (this can be neglected for quasars if we consider 2: > 1), 
their random velocities and redshift errors. This effect is es- 
pecially noticeable for quasar pairs with the projected linear 
separations < 2 Mpc, where the SPS S data are expected to 
be incomplete l|Hennawi et al.ll2006l ). On larger scales the 
2pCF profile is flatt ened along the line of sight ('The Bull's 
eye' or'Kaiser ( f 987) effect) due to the gravitational infall to 
density inhomogeneities; this kind of redshift-space distor- 
tion dominates on the linear scales. These effects are param- 
eterized by line-of-sight pairwise velocity dispersion {w^Y^^ 
and infall parameter /3. It is worth to note that the non- 
Hubble motions and redshift errors are not the only sources 
of redshift-space distortions. One more effect of geometric 
flattening that can lead to distortion of 2pCF can be due 
to a wrong choice of cosmological parameters Q.m, ^a; this 
provides an additional tool for estimat ion of these parame- 
ters b y means of a geometrical test of lAlcock fc Paczvnskil 
l|l979[) . The redshift-space 2pCF of quasars can be in prin- 
ciple used for estimation of all parameters (/?, (w^)^'''^ and 
cosmological ones) simultaneously. But due to a degener- 
acy between the geometric distortions and the redshift-space 
distortions ( see, e.g., iHoyle et al.l (l2002i)') this problem is 
complicated. iHovle et all l|2002l )"and lda Angela et al.l (120051 ) 
proposed a method that all ows to break this degenerac y by 
means of combination of the lAlcock fc Paczvnskil l| 19791) test 
with that based on evolution of the quasar clustering ampli- 
tude, which has a different dependence on /3(z) and Qm{0). 

As we cannot estimate the non-Hubble motion of each 
quasar independently, the effects of the redshift-space dis- 
tortion are at present the only source of statistical informa- 
tion about proper velocities of quasars. On the other hand, 
these effects prevent direct determination of 3D 2pCF that 
involves line-of-sight distances between objects determined 
from z- measurements. This urges us to use only the pro- 
jected distances a (i.e. orthogonal to the line of sight, which 
may be considered as independent of proper motions) to 
determine the projected 2pCF and then to restore the real- 
space 2pCF. It is well known that such reconstruction of the 
real-space 2pCF from projected one is mathematically ill- 
posed problem. However, most authors avoid this difficulty 
using a concrete functional form of 2pCF. Typically 2pCF 
is represented in a power-law form ^(r) = (ro/r)^, though 
it is clear that 2pCF slope and correlation len gth may be 
differe nt on different interval. For example, ida Angela et al.l 
showed that double power-low model gives a good 
fit to 2pCF of quasars distribution with ^(r) = (e.G/r)^ **^ 
over scales of 1 < r < 10 /i"^ Mpc and ^ = (7.25/r)^'^° 
over scales of 10 < r < 40 Mpc for a quasar sam- 



ple from 2QZ catalogue with redshifts within the range 
0.3 < z < 2.2. Using quasar catalogue of spectroscopi - 
cally confirmed quasars (|Schneider . Hall fc Richards' 200'n) 
based on the 5th data release of SDSS Ross ct al. (2009|) 
obtained ro = 5.45to j5 h''^ Mpc and 7 = 1.90^0.03 over 
scales of 1 < cr < 130 Mpc for low-redshift (0.3 < 
z < 2.2) quasars, and for high-re dshift (2.9 < z < 5.4) 
quasars from the same catalogue IShen et al.l (|2007l ) ob- 
tained ro = 15.2 ± 2.7 Mpc and 7 = 2.0 ± 0.3 
over scales 4 < cr < 150 Mpc. Note that similar 
large corr elation length was o btained for radio-sources by 
iBlakc, Mauch fc Sadleij (l2003l) and for X-ray sources by 
Ipiionis et all l|2004l ) ajid iBasilakos et all (|2004l ): at least a 
part of both of these object types are considered to be 
AGNs. In addition to r edshift depen d ence o f the quasar 
clustering confirmed by Groom et al. ' J200|ii_ iMyers et slI] 
(2006, 2007);_ Porciani, Magliocchct ti fc Norber g (200^ 
IPorciani fc Norberel (l2006i): Ida Angela et al.l (I2008D some 



authors dPorciani fc Norberd I2OO6I : Ida Angela et all I2OO8I : 
iMountrichas et aL 20091 ) argue that quasar clustering de- 
pends weakly on luminosity. Hence speaking about repre- 
sentation of 2pGF in power-law form in the present paper 
we always mean that 2pGF is averaged over quasars with dif- 
ferent luminosities and over certain redshift range and dis- 
tance intervals described below. In particular, when dealing 
with the bias estimates we confine ourselves to the intervals 
where our results may be compared with the linear theory 
of matter inhomogeneity. 

The most complete full investigation on non-Hubble 
motions of quasars together with cosmological parame- 
ters was done by 2 dF team for 2dF QSO Survey (see 
iGroom et all (|2005|): [ da Angela et al.l (|2005|): iHovle et all 



( 20021) : lOutram et al.l (|2004l )). E.g. Ida Angela et al.l (|2005D 

obtained = 0.35lg;l| and /3 = 0.50lg:l| (for the sam- 
ple mean redshift z = 1.4). In these studies the value of 
pairwise velocity dispersion of quasars was usually treated 
as fixed parameter and its value was chosen either from 
Hubble volume simulations or following previous results on 
galaxy pairwise velocity dispersion assuming that this value 
does not evolve strongly with time. For galaxies, we have 
more definite results for the velocity dispersion obtained 
from the redshift-spac e distortions. The first estimate of 
I Davis fc PeeblesI (|l983l ). 340 ± 40 km/s, remained a canoni- 
cal one almost a decade. However, later IMo. Jing fc Bornerl 
( 1993 ) showed that the value of pairwise peculiar veloc- 
ity dispersion depends strongly on the sample. Later this 
res ult was confi r med by values 540 ± 180 km/s, given 
by iMarzke ei~sl\ (| 19951 ) for the Second Genter for Astro- 
physics Galaxy Redshift Survey (CfA2) and the South- 
ern Sky Redshift Survey (SS RS2), and 416 ± 36 km/s, 
given by iRatcliffe et al.l (| 19981 ) for Durham/United King- 
dom Schmidt Tel escope (UKSTj Galaxy Redshift Survey 
and others. Ratcl iffe et al.l (|l998h also presented the value 
P — 0.48 =b 0.11, which agree s well with later results, such as 
the results of lHawkins et all (|2003l ) for 2dF Galaxy Redshift 
Survey. A similar value, 330 km/s, was obtained for lumi- 
nous red galaxies (LRG ) from the 2dF -SDSS LRG and QSO 
(2SLAQ) survey bv .Ross et al.l (|2007| ). 

In this paper, we present new results on the 2pCF slope 
and correlation length, infall parameter /3 using the last 7-th 
Data Release of SDSS. In Section [2] we describe the quasar 
samples and technique for random catalogue generation we 
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Figure 1. Sky coverage of our full quasar sample in equatorial 
coordinates. 

used in our study. In Sections[3]and[4]we present the results 
on parameters of the real-space and redshift-space 2pCFs 
and techniques for their estimation. Section |4] also contains 
estimations of the infall parameter /3 and quasar bias b from 
correlation function analysis on linear scales larger than 
10 Mpc. In Section [5] we discuss possibilities to use 'Fin- 
ger of God' effect to estimate the random quasars velocities, 
which does not correlate with density inhomogeneities. Be- 
cause determination of random pairwise velocity dispersion 
{w'^)^^'^ of quasars needs knowledge of 2pCF in the region 
r < 2 Mpc, which is not well studied, we considered differ- 
ent simpUfied models of 2pCF to study their outcome for 
{w^y^^. Finally, in Section [6] we sum up the results. 

Throughout the paper we use the spatially flat cosmo- 
logical mo del with the combined estimates f^A = 0.726 and 
h = 0.705 (|Komatsu et a l. 2009), which agrees with recent 
estimate from seven-y ear Wilkinson Mic rowave Anisotropy 
Probe (WMAP) data (|Larson et al.ll2010l ). 



2 DATA 

2.1 The sample 

Our sample is taken from the 7th Data Rele ase of the 
Sloan Digital Sky Survey (|Abazadiian et al]|2009l ). contain- 
ing about 100,000 quasars. The redshift range used in our 
analysis is 0.8 ^ z 2.2 and as the sky coverage of SDSS 
contains a one big piece and three narrow near-equatorial 
'stripes', we excluded these stripes to reduce boundary ef- 
fects. As the data was taken from the list of all objects clas- 
sified as quasars, which is not a vetted catalogue, we firstly 
excluded 'bad' objects: objects with failed astrometry and 
photometry data. The second step was to exclude objects 
with unreliable redshifts. The redshifts in SDSS are mea- 
sured mainly using two techniques, emission-line measure- 
ments and cross-correlation, either both or one of these. We 
excluded objects with the redshift status 'INCONSISTENT' 
and 'FAILED'. The first status means that the redshift was 
measured by two techniques but the results were inconsis- 
tent. The second one means that the redshift measurement 
failed. The objects with redshift error > 0.01 were also ex- 
cluded. 



Figure 2. Redshift distribution of our quasar samples: full 
(solid line), uniform (dashed), low-reddening (dash-dot), high- 
reddening (dash-dot-dot-dot), good (long dashes) 



-50 r 

-29 - 



-28 - 




-23 L , , I , , , I , , , I , , , I , , , I , , , I , , ^ 

0.8 1,0 1-2 1,4 1,6 1-8 2-0 2-2 

redshift, z 

Figure 3. Luminosity — redshift distribution for uniform sample. 

The third step was examination of all objects left af- 
ter the previous steps by eye (photometry and spectroscopy 
data) using the SDSS data base. During these examination 
466 objects were excluded as the ones having too faint and 
noisy spectra, or being fault 'double' objects (i.e. single ob- 
jects processed twice and thus appeared twice in the list 
with similar coordinates), or being artifacts of observations, 
or being stars (we have found 5 stars in this sample). 

Our final sample, which we call full, contains 52303 ob- 
jects and has the mean redshift z = 1.47. The sky cover- 
age of the sample and its redshift distribution are shown 
in Figs. [1] and [2] (solid line). The redshift-luminosity distri- 
bution of this sample is shown in Fig. [3l where the effect 
of nii — 19.1 magnitude limit can be seen. Here the values 
of absolute magni tudes in i-band, Mj, a re calculated using 
./f-correction from [Richards et al.l l|2006l ). 

2.2 Sample inhomogeneity 

As one can see from Fig. [T] the sky coverage by our 
sample in not heterogeneous: some parts seem to be signif- 
icantly denser then others. This is due the fact that these 
parts were twice spectroscopically covered and can influence 
our results. First of all we note that the algorithms we use 



4 G. Ivashchenko, V. I. Zhdanov, A. V. Tugay 



to estimate the 2pCF parameters (sections 3.1 and 4.1) are 
just worked out in order to r educe the selection effects due 
to unphysical inhomogcneitv |lvashchenko fc Zhdanovl2010l : 
IZhdanov fcl^ashche nko 20m. 

Nevertheless, to test the influence of inhomogeneity we 
constructed 50 homogeneous samples in the following way. 
The whole sky area covered by our sample was divided into 
2.5° X 2.5° patches and the density of each patch, mean 
density (n — 7.25 objects per square degree) and its rms 

= 2.55 objects per square degree) were calculated. All 
the objects from the patches with density exceeding the 
mean one more than 0.5rms were excluded (see Fig.|4|. Then 
we generated 50 new homogeneous samples by putting the 
excluded objects in their places, but choosing (in a random 
way) only those making the density of the given patch equal 
to the mean density. The sky coveredge of one of these sam- 
ples is presented in Fig. [5] Each new test sample has 35643 
objects and a density n = 6.35 ± 1.19 objects per deg^. Ap- 
parent decrease of the density with increasing declination is 
the result of projection of the celestial sphere onto the plane. 



2.3 Tilable targets 

Most quasar candidates in the SDSS are selected 
based on their locations in multidimentional SDSS colour 
space and their cross-identification with radio sources 
from the Faint Im ages of the Radio S ky at Twenty-cm 
(FIRST) survey (|Richards et al.l 120021 '). Supplementing 
this primary quasar sample are quasars targeted by the 
GALAXY, X-RAY, STAR and SERENDIPITY SDSS 
software packages; no attempt at completeness was made 
for the last three categories. Thus, the sample we use 
is not complete and the possible effect of this incom- 
pleteness on our results has to be investigated. For this 
purpose we constructed the sample of the so-called 'tilable' 
objects (uniform) from our sample. According to SDSS 
glossary web-page, tilable targets are supposed to have 
as closed to uniform completeness as possible. These are 
targets with primTarget flag 'QSO_HIZ', 'QSO.CAP', 
'QSO_SKIRT', 'QSO_FlRST_CAP', 'QSO_FlRST_SKIRT', 
'GALAXY_RED', 'GALAXY', 'GALAXY_BIG', 
'GALAXY_BRIGHT_CORE', 'STAR_BROWN_DWARF', 
and second target flag 'HOT_STD'. The number of objects 
in the subsample constructed following this criterion is 
37290. One can see that the sky coverage by the uniform 
sample presented in Fig. [6] is close to uniform and close to 
the sky coverage by our homogeneous samples (Fig. [5]). The 
redshift distribution of this uniform sample is presented in 
Fig. [2] (dashed line). 
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Figure 4. Sky coverage of the full sample with excluded patches 
with density higher than O.Srms. 
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Figure 5. Sky coverage of one of the homogeneous samples. 

0.2603) one. Values of E{B — V) were obtained from the 
reddening in each SDSS band given in the catalogue and us- 
ing the relative extinction for different bandpasses from Ta- 
ble 6 in Schlegel. Finkbeiner fc David (|l998l ). The numbers 
of objects in subsamples are 30188 and 22115 for low and 
high reddening correspondingly. The redshift distributions 
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2.4 Reddening 

As the main criterion of quasar candidates selection in 
SDSS is based on their magnitudes and colors, that have 
been corrected for Galactic e xtinc tion using the maps of 
ISchlegel. Finkbeiner fc David (|l998l ). any systematic errors 
in the reddening model can induce a dditional effects on clus- 
tering results (see IRoss et al] (|2009|^ for discu s sion o f possi- 



ble effects). That is why following Ross et al] (|2009l ) we de- 
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vided our sample into low reddening (0.0028 < E{B -V) ^ pi g_ gj^ coverage of the umform sample. 

0.0217) sample and high reddening (0.0127 < E{B - V) ^ 
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of these subsamples are presented in Fig. [2] (dash-dot for low- 
reddening sample and dash-dot-dot-dot for high-reddening 
one). 

2.5 Bad fields 

Since the quasar selection algorithm involves colours 
and magnitude all photometric errors caused by e.g. pho- 
tometric callibration can induce errors to this selection and 
result in additional inhomogeneity of the sample. For veri- 
fication of these effect we constructed a good sample using 
information about the data quality for each field (an image 
in all five bands with approximate dimentions of 13' x 10') 
from a field quality table, called RunQA. To construct the 
good sample we used the FieldQAU number which describes 
the overall quality of a field and rejected the objects with 
its value equal to 'bad'. This overall quality determination, 
FieldQAU, is based on the (dereddened) principal colors of 
all the stars in the field vs. the Galactic value of this quan- 
tity, the PSF quality, the difference between Aperture and 
PSF magnitudes for the same stars in the field, and the see- 
ing. We also excluded the objects with rerun number equal 
to zero, as there are no fields with such rerun number in 
FieldQAU table. The number of objects in the good sample 
is 50420. Its redshift distribution is presented in Fig. [2] (long 
dashes) . 

2.6 Fibre collisions 

One another peculiarity of our sample is related to the 
so called 'fibre collisions' effect, which is inherent for mul- 
tifibre spectrographs. It means that two spectroscopic tar- 
gets on the same plate cannot be measured simultaneously if 
they are closer than some distance (55" for SDSS); it is only 
possible to obtain spectra for both such objects if they are 
covered by more than one plate (Schneider, Hall & Richardi 
l2007h . Angular separation 55" corresponds to 0.7— 1.5 Mpc 
projected distance for our redshift interval 0.8 — 2.2. That is 
why we avoided this distance range when calculating the 
parameters of 2pCF. Of course this effect affects mostly 
the precision of our random velocities dispersion estimations 
due to the small sample of pairs with projected separation 
< 2 M pc. But in this case we cannot follow e. g. IRoss et aU 
l|2009t ) supplementing their sample with quasar candidates 
from SDSS Quasar photometric catalogue with photometric 
redshifts ((Richards ct al. 2 009i) due to the fact that the er- 
rors of photometric redshifts are about Az = 0.3(2(j) which 
is much larger than the restrictions we put on our sample. 

2.7 Random samples 

An important point in reconstruction of the correla- 
tion function of extragalactic objects is generation of the 
random sample, which should reproduce a random distri- 
bution of objects with the same redshifts and angular dis- 
tributions inherent to the initial sample as much as possi- 
ble. For this pu r pose an imaging mask may be used (see 
e.g. iMvers et al.l (|2006l )'). In the present work we used the 
p er mutatio n met h od of random catalogue g e neration (see 
e.g. lOsmerl (|l98ll ): IZhdanov fc Surdeil (|200ll ): ICroom et all 
(|2005h ') that does not require any technical information 



about the survey. This method uses the random permuta- 
tion of objects' redshifts from the initial list of quasars while 
keeping their coordinates unchanged. As a result one obtains 
a randomized list of objects with the same redshift distribu- 
tion and the same distribution over the sky as in the initial 
sample. At the same time this randomization destroys cor- 
relations between objects of the initial sample: those objects 
which were clustered in the initial sample appear to be not 
clustered in the randomized sample as they have most prob- 
ably different redshifts and thus are located far from each 
other. 



3 REAL-SPACE CORRELATION FUNCTION 
3.1 The technique 

As peculiar velocities of objects distort the distance be- 
tween objects along the line of sight, the parameters of the 
real-space 2pCF of quasars can be reconstructed from the 
projected 2pCF. We use this standard approach, but tak- 
ing into account some specific features of the permutation 
method of random catalogue generation. 

2pCF ^(r) determines the probability dPi to find a 
neighbour of the object Qi from the sample at co moving 
distance r in a volume element dV as (|Peebleslll98"ol ) 

dn = n,[l + ({r)]dV, (2) 

where Ui is a number density of i-th object in the neighbour- 
hood of Qi. Here we must specify the definition of the dis- 
tance between neighbours, which can be defined differently 
in curved space-time. For the distance scales under consid- 
eration this may be done with a sufficient accuracy using 
local Lorentz frames (|Misner. Thorne"fc" Wheeleij[l973l ). For 
definiteness we use the Fermi reference frames and normal 
Fermi coordinates (Synge 1960; Misner et al. 1973); this is 
especially convenient because in case of spatially flat cosmo- 
logical background this leads to distance relations analogous 
to that of Euclidean space (see Appendix) . 

We suppose that the real-space 2pCF of quasars has 
power-low form ^(r) = (ro/r)^ for distances r < L and 
r > Lo, where L is some characteristic comoving scale where 
2pCF approaches to zero. We also make the following as- 
sumptions. First, we neglect variations of Ui on the scales 
< L. In fact, we also considered models that take into ac- 
count variations of Ui with redshift on this scale and found 
the resulting difference to be not statistically significant, as 
it can be expected from very beginning. Second, we assume, 
as usual, that quasars distribution is locally isotropic at any 
point of the real-space, i.e. it depends only on the distance 
between the objects. Third, we neglect variations of redshift 
distribution of the sample in different sky regions. 

Further we pass to projected (comoving) separations. 
The only thing th at makes our consideration different from 
the standard one (|Peebleslll98"oh is related to our method of 
the random catalogue generation. When counting the num- 
ber of pairs with projected separations within the range 
[a, a -\- da] , we include only pairs with l ine-of-sight separa- 
tions H that satisfy the condition H < L (|Zhdanov fc Surdeil 
I2OOII : [Zhd anov fc lvashchenk3 120081 ). (Here we denote the 
distance along the line of sight as H to distinguish it from 
the number tt.) This number of pairs should be compared 
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to analogous number in the randomized sample. If, e.g., we 
do not restrict 11 and include all pairs with the above pro- 
jected separation, there obviously would be no difference in 
these numbers between the initial sample and the random- 
ized one obtained by permutation of redshifts. On the other 
hand, we expect that almost all the pairs that have a notice- 
able correlation are concentrated within the region 11 < L. 
Note that though at this point we use the observational val- 
ues of n that involve contribution of the proper velocities, 
corresponding correction due to velocities here is negligible, 
because it affects only larger distances ~ L. 

The total number of neighbouring quasars with the 
projected distances from a given quasar within the range 
[a, a + da] and line-of-sight separations 11 < L is obtained 
by summing over all quasars from the sample 



DD{a, da) ^ ^ 



rii 



2L + 



n2 du 



dS, (3) 



where the sum is performed over all objects in the sample, 
dS = 2TTada. This number is twice larger than the number of 
pairs with separation from the above range, if the quasars 
and the neighbours are chosen from the same sample. In 
order to study boundary effects we also considered the case 
when the quasars are chosen from a subsample obtained 
from the initial sample with the exception of objects near 
the boundary, and the neighbours being looked for in the 
whole initial sample. 

In case of the randomized redshifts the corresponding 
number of neighbours is 

RR{a, da) = 2L^n^dS', 

i 

where n'^ is the number density in the neighbourhood of 
i-th object in the randomized sample. As the redshift dis- 
tribution of objects in initial and random samples are the 
same, than one estimates "^rn ^ '^n'^ with a relative error 



~ 1/VN, where A'' is the number of objects in the sample. 
Therefore 



DD{a, da) 
RR{a, da) 



1 + 



2L 



n2 dn, 



(4) 



where we extended the integration limits to infinity. 
Note that possible difference of the ratio X] '^«/ 

i i 

from unity may be due to variations of z-distribution of 
quasars in the different sky regions of the real survey. Cor- 
responding correction was analyzed and found to be statis- 
tically insignificant (jZMaiiov & Ivashchenko 20o3). 

In fact we used not exactly Eq. Q but correspond- 
ing relations for separations in finite intervals [a, a + Aa] ; 
typically we used Aa = 1 Mpc. Then using the power- 
low form of 2pCF we have the follo wing relations for fit - 
ting parameters 7 and rp (cf., e.g.. [ Shanks et al.l (|l980h : 
ICollins. Hevdon-Dumbleton fc MacGillivravl jigsgh ): 



RR{a, Aa) 
where 



B 



{a + Aa 



,3-7 



3-7 



(3-7)(a + Aa/2) 



(5) 



B = 



£7, Mpc 



Figure 7. Dependence of normalized neighbour numbers dd (for 
the initial sample) and rr (for the random sample) on projected 
distance a. Top; full sample, bottom: test, described in text. Cir- 
cles denote initial sample, triangles denote random sample. 



The errors of all the parameters in the present study 
have been es timated using the 'jackknife' method (see, e.g., 
iMilleil ||1974 ) for review), which has been successfully ap- 
plie d to clustering prob lems in extragalactic astronomy (see, 
e.g.. lMvers et al.) (|2006[ )). We divided our initial sample into 
Na ~ 50— 100 equal stripes along declination; owing to large 
number of objects in each stripe, the stripes may be viewed 
as practically independent ones (here we neglect the bound- 
ary effects). Then new Ns subsamples have been formed: i-th 
subsample is unification of all the stripes except i-th stripe. 
Jackknife dispersion estimate of any parameter x calculated 
in this treatment is then 



7V- 



iV 



1 ^ 



(6) 



2Lr(; 



where Xi is the corresponding value of this parameter for 
i-th subsample and x is either the value of x for the initial 
sample or it is average over Ns subsamples. 

3.2 Results 

We determine the distances between quasars using the 
locally Lorentz frame of the first qua sar (see Appendix ). 
This method differs, e .g., from that bv iHovle etall (j2002l ): 
Ida Angela et al.l (|2005l ). 

Our sample covers only a part of the sky; thus the 
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Table 1. The parameters of the real-space 2pCF for different cases described in text and the results of other authors for comparison. 
Here DR5Q stands for quasar catalogue by Schneider et al. (2007). 



dist. range, h ^Mpc z range mean z rg, h ^ Mpc 7 jd.o.f. sample authors 















5.85 ± 0.33 


1.87 ±0.07 


1.054 


full 
















5.83 ± 0.36 


1.87 ±0.08 


2.084 


test 








df 








5.83 ±0.36 
5.76 ±0.84 


1.87 ±0.08 
1.78 ±0.26 


2.084 


test 
homogeneous 




fT 


= 1 - 


35 


0.8 


- 2.2 


1.47 


6.08 ± 0.59 
5.84 ± 0.39 
5.63 ± 0.55 
5.73 ± 0.57 


1.88 ±0.11 

1.87 ±0.09 

1.88 ± 0.12 
1.85 ± 0.11 


0.789 
1.231 
0.846 
1.025 


uniform 
good 
low reddening 
high reddening 


present paper 


r 


= 1 - 


20 


0.8 


~ 2.1 


1.47 


4 s+O-9 
^•^-1.5 


1.5 ±0.2 




2QZ 


Porciani et al. (2004) 


r 


= 1 - 


10 


0.3 


- 2.2 


1.4 


6 0+° '^ 
"■"-0.7 


1.45 ± 0.27 




2QZ 


da Angela et al. (2005) 


r 


= 10 - 


- 40 


0.3 


- 2.2 


1.4 


7.25 


'^■''"-0.03 




2QZ 


da Angela et al. (2005) 


a 


= 4 - 


150 


2.9 


- 5.4 




15.2 ± 2.7 


2.0 ±0.3 




SDSS DR5Q 


Shen et al. f2007) 


a 


= 1 - 


130 


0.3 


- 2.2 


1.27 


5 45+0.35 


1 90+°" 




SDSS DR5Q 


Ross et al. ("2009) 



boundary effects can influence the result because the num- 
ber of neighbours for objects lying close to the boundary 
is less then for those lying far from it. To estimate this 
influence we made calculations in two different ways. In 
the flrst one we calculated the number of neighbours for 
all the quasars in the sample. In the second one we took 
into account the neighbours of the objects lying at least 
5 degrees far from the boundary (we called this test sam- 
ple). Fig. [7] shows the surface number density of neighbours 
dd — DD/[nAa{2a + Aa)] in the initial sample (circles) and 
rr = RR/[TTAa{2a + Aa)] in the random one (triangles) cor- 
respondingly, where we choose the bin size Act = 1. Values 
RR are the mean numbers over 50 realizations of random 
sample, error bars represent Ict jack-knife errors. 




As can be seen from Fig. [T] in practice the surface den- 
sity of neighbours for the random catalogue does not depend 
on the distance (except for the distance range < 1 — 2 Mpc), 
which is expected for a random distribution of objects. The 
gap at ~ 1 — 2 Mpc can be explained by a deflcit of close 
pairs in the initial catalogue. Note that angular distribution 
of objects in our random samples is the same as in the initial 
one (see Sec. 12. 7p . This deflcit is a result of the fibre colli- 
sion effect (Sec. 12. 6p . One can see that the number densities 
for initial and random samples coinside up to 40 — 50 Mpc, 
thus we use L = 50 Mpc (or 35 Mpc). And we also put 
Lo = Ih-^ Mpc because of the fiber collision effect. The 
quasar projected 2pCF Wp{a) divided on ct in logarithmic 
scale for the full sample is presented in Fig. |S] with the the 
best fit single power-low over the range 1 < a < 35 Mpc. 

The same calculations were provided for 50 test samples 
(see Sec. 12. 2p . The scattering of the resulting parameters is 
presented in Fig. [9] with the mean values presented in the 
second row of Table 1, where errors are rms. Results for 
homogeneous, uniform, good, low reddening and high red- 
dening samples are presented in rows 3 — 7 of this table. For 
these samples the errors are jackknifes. Comparison of the 
projected 2pCF for the initial sample with uniform, good, 
low reddening and high reddening samples are presented in 
Figs. [T0TT2I 



Figure 8. The quasar projected 2pCF, Wp, for full sample. The 
quoted errorbars are jackknifes. The solid line showes the best fit 
single power-low with ro = 5.85 /i""*^ Mpc and 7 = 1.87 over our 
full range 1 < cr < 35 Mpc. 
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Figure 9. Correlation lengthes and slopes of the real-space 2pCF 
for 50 homogeneous samples. 
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Figure 10. The quasar projected 2pCF, Wp, for initial (squares 
and solid line) and uniform (triangles and dashed lines) samples. 
The quoted error bars are jackknifes. 




3.3 Discussion 

As one can see from Figs. 11011121 and Table 1, the pro- 
jected 2pCFs and parameters of the real-space 2pCF for all 
samples agree well with each other within the errors that 
are slightly larger for the samples with smaller number of 
objects. Hence, firstly, we can neglect boundary effects for 
our sample, that is expected due to the sample size being 
several times larger than the distance range we consider. 
Secondly, we can also neglect all the inhomogeneities of the 
sample caused by the non-uniformity of the sky coverage as 
a result of the spectroscopic survey, reddening and the con- 
ditions of observations. The last five rows of Table 1 show 
comparisons of our results with the results of other authors. 
We can see that our results agree well w ith early results 
for the 2QZ catalogue by Porciani. Mag liocchetti fc Norberd 
(|2004l ) within Icr for rp and 2a for 7. The best agreement 
is with the results by IRoss et al.l (|2009l ) for the 4th edition 
of SDSS quasar c atalogue (compiled from the 5th d ata re- 
lease of SDSS by ISchneider. Hall fc Richard j (|2007l ')'). but 
note that they used the wider dis tance range. The 2pCF 
parameters from IShen et al.l l|2007l ) are included to the ta- 
ble only for completeness of known results on the quasars 
2pCF in the real-space as they correspond to larger red- 
shi fts. Our results c annot be directly compared with that 
bv lda Angela et all (|2005l ) for 2QZ catalogue, though there 
is rather a good agreement if we average their values over 
both distance ranges 1 — 10 and 10 — 40 Mpc. Note that 
our redshift range is smaller, but the fraction of quasars 
contained within the 0.3 < z < 0.8 range is small compared 
with those in the 0.8 < 2 < 2.2 range. Finally, we note 
that comparison of 2pCF parameters on diflerent distance 
intervals shows that the power-law form of 2pCF must be 
replaced by more complicated one for larger distances. 



Figure 11. The quasar projected 2pCF, Wp, for initial (squares 
and solid line) and good (triangles and dashed lines) samples. The 
quoted errorbars are jackknifes. 




Figure 12. The quasar projected 2pCF, Wp, for initial (squares 
and solid line), low-reddening (triangles and dashed lines) and 
high-reddening (diamonds and dotted lines) samples. The quoted 
errorbars are jackknifes. 



4 REDSHIFT-SPACE CORRELATION 

FUNCTION AND INFALL PARAMETER 

4.1 The technique 

2pCF in the redshift-space is not spherically-symmetric 
because of the eflects of velocities mentioned in Sec. [T] Here 
we pay attention mainly to the scales larger than 10 Mpc, 
where we expect the density contrast to be sufficiently small. 
In this region the effect of bulk gravitational infall leads to 
flattening of redshift-space 2pCF contours and to certain 
enhancement of the monopole part ^o{s) of this 2pCF, which 
is obtained by averaging over all angles. On the other hand, 
in this region we neglect 'Finger of God' effect due to random 
velocities. 

The relation of the real-space 2pCF to the redshift- 
space 2pCF is well known within the l inear theory of 
cosmological perturbatio ns (|Kaiseil Il987l : iHamiltonl Il992l : 
iMatsubara fc Sutdll99d ). The monopole (i. e. angle aver- 
aged) part of 2pCF is proportional to the real-space 2pCF 
iKaiseil [19871 ') and therefore it also has a power-law form 
Co(s) = {so/sy" with the same slope 7^ = 7, but with dif- 
ferent correlation length: 



So = ro 



, 2/3 ^2 
1 -I- — + — 
3 5 



l/7a 



(7) 



This formula is well-grounded in the region r > 10 Mpc, 
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where the hnear theory is apphcable. Here the parameter 
P = f{^M,z)/b appears as some consequence of the large- 
scale infall onto matter overdensities, h is the bias parameter 
in Eq. ([l]). The function f{Q.M,z) is express ed by means o f 
the growing mode of the density contrast JPeeblesI 198(J). 
There is a simple analytical approximation l|Carroll fc PressI 
1 19921 ) for this function in the case of the spatially-flat ACDM 
cosmological model: 



M,Z) 



nM(l + 2)3-Hl-f^i\ 



4/7 



(8) 



Therefore we may proceed in two ways: either (i) by 
taking the slope 7s = 7 from the results of Sec. [3] dealing 
with fitting of so only, or (ii) by independent determination 
of 7s, So. 

Like in Sec. [3] one can determine the total number of 
neighbours of the i-th quasar from the whole sample in a 
spherical layer [s, s + As] as 

t-As 



DZ)(s, As) = 



47r 



1 + ^ / as')s'^ds' 



where AV = 47rA(s^ -I- sAs -I- (As)^/3). The similar estima- 
tion for random catalogue, which is considered to represent 
random spatial distribution of objects with no clustering, is 

RR{s,As) = iTv^n'iAV. 

i 

Assuming X] ~ X] ^^'^ taking into account the power- 

i i 

low form of the monopole part of 2pCF 5o(s) — {so/sY'" , we 
have the relation for fitting: 



DD(s, As) 



(s + As) 



,3-7, 



RR{s,As) 3-7s ' As(s2-fsAs-HAs2/3)' 

In order to determine the infall parameter /3 on the 
way (i) we use Eq. ((9]) to find so, with the slope 7^ = 7 
being known from the results of the Sec. O Then we derive 
P by solving the Eq. (O with known so, ro, and then we 
estimate bias b for our redshift interval using Eq. (|8} . These 
operations deal with the region s > 10 Mpc. In addition, 
mainly for comparison with the results of the other authors, 
we also considered way (ii) with independent determination 
of 7 and So for different separation intervals. 

Note that the method we use for estimation of the in- 
fall para meter is similar to the so-called J3(s)/J3(r)-method 
(see e.g. iRatcliffe et al.l (|l998l )). which uses the volume in- 

X 

tegrals J3{x) = J ^{y)y^dy. Using of these integrals allows 



to smooth the influence of small sizes of the samples. Our 
approach has the same advantage as it includes the 2pCF 
parameters which are the averaged characteristics of 2pCFs 
within the whole distance range. 



4.2 Results and discussion 

We estimated the 2pCF parameters within the whole 
range s = 1 — 35 Mpc and separately for ranges 

1-10 Mpc and 10 - 35 h"^ Mpc. The resuhs are pre- 
sented in Table 2 and Fig. 13. The errors are calculated using 
the jackknife technique. We also found the same parameters 
for dist ance range 1 — 25 ft" ^ Mp c for compa r ison w ith the re- 
sults of lCroom et~aLl (120051 ) and lRoss etal] |200i) presented 




Mpc 



Figure 13. The quasar redshift-space 2pCF, for full sample. 
The quoted errorbars are jackknifes. The best fit single power-low 
over our full range 1 — 35 Mpc and over ranges l — 10h~^ Mpc 
and 10 — 35 Mpc with parameters presented in Table 2 are 
shown with solid, dotted and dashed lines correspondingly. 
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Figure 14. Correlation lengthes and slopes of the redshift-space 
2pCF for 50 homogeneous samples. 



in the last 3 rows of this table. As one can see our parameters 
So and 7s within the d istance range 1 — 25 Mpc agree 
well with the results of lCroom et al. I HooF) andU OSS et al.l 
(2009) for 2QZ ca talogue and the 4th edition of th e SDSS 
quasar catalogue (|Schneider. Hall fc Richardsl |2007| ) corre- 
spondingly. An d for distance range 1 — 10 Mpc our re- 
sults agree with lCroom et all (|2005l ) within 2a for the slope 
and la for the correlation length. 

Comparing the parameters for different distance ranges 
we can conclude that on the large scales the redshift-space 
2pCF is steeper than on the small scales while the correla- 
tion lengthes do not differ significantly within la. This dif- 
ference in the slope of the redshift-space 2pCF of quasars o n 
different scales was pointed out earlier (|Croom etal .1120051 '). 
It is important to note that the results for the slope on 
10 — S5h~^ Mpc interval fairly well agree with the results 
for 7 in the real space, as it must be. 

As discussed in Section 4.1, we have used the correla- 
tion length ro and the slope 7 of the real-space 2pCF ob- 
tained in Section |3l We have then obtained the redshift- 
space correlation length so = 7.27 ± 0.37 Mpc for the 
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Table 2. The parameters of the redshift-space 2pCF and the results of other authors for comparision. Here DR5Q stands for quasar 
catalogue by Schneider et al. (2007). 



dist. range, 


fe-^Mpc 


z range 


mean z 


So, h~i Mpc 


7s 


xVd.o.f. 


sample 


authors 


s = l 


- 35 


0.8 - 


2.2 


1.47 


5.83 ± 0.47 
5.95 ±0.86 
5.33 ±0.51 
5.76 ±0.54 
5.45 ±0.96 
4.56 ±0.55 


1.49 ± 0.12 
1.44 ±0.21 
1.35 ±0.08 
1.46 ±0.12 
1.48 ±0.24 
1.24 ±0.11 


1.440 

1.473 
1.363 
0.779 
1.892 


initial 
homogeneous 
uniform 
good 
low reddening 
high reddening 


present paper 


s = 1 
s = 1 
s = 10 


- 25 

- 10 

- 35 


0.8 - 


2.2 


1.47 


5.97 ±0.51 
6.43 ± 0.63 
7.37 ±0.81 


1.37 ±0.13 
1.21 ±0.24 
1.90 ±0.24 


1.060 
0.456 
1.636 


initial 


present paper 


s = 1 
s = 1 
s = 1 


- 25 

- 10 

- 25 


0.3 - 
0.3 - 
0.3 - 


2.2 
2.2 
2.2 


1.35 
1.35 
1.27 


5.48_(, 4g 
o Ho+0.43 
•^•88-0.53 
5.95 ±0.45 


1.20 ±0.10 




2QZ 
2QZ 
SDSS DR5Q 


Groom et al. f2005"l 
Groom et al. ("20051 
Ross et al. (20091 



Table 3. The values of infall parameter /3 and bias parameter bg and the results of other authors for comparison. Here DR5Q 
stands for quasar catalogue by Schneider et al. (2007), DR4phot stands for a sample of photo metrical l y clas sified quasars fro m the 
4th D ata Re lease of SDSS . The numbers in brackets correspond to methods pr oposed by Hovl e et al.l | |2002|1 - [l],|da A ngela et al.l 
1I2OO5I I - f2], lGroom et al.l ll2005h - [3], and lPorciani. Magliocchetti fc Norberd ||2004I 1 - [4]; see explanations in the text. 



mean z bq f} method sample authors 



0.55 


1.4 ±0.2 


0.55 ± 0.10 


mean of [1],[2] 


2SLAQ+2QZ+SDSS DR5 


Mountrichas et al. (20091 


1.27 


2.06 ± 0.03 


0.43 


[3] 


SDSS DR5Q 


Ross et al. (20091 


1.35 


2.02 ± 0.07 


0.44 ± 0.02 


[3] 


2QZ 


Croom et al. (20051 


1.40 


1.50 ±0.20 


o.6o«:lt 


[2] 


2QZ+2SLAQ 


da Angela et al. (2008) 


1.40 




n 09+0.09 
'^•'^^-0.11 


i{s)/i{r) 


2QZ 


da Angela et al. (2005) 


1.40 






[2] 


2QZ 


da Angela et al. (20051 


1.40 


2.41 ±0.08 




[3] 


SDSS DR4phot 


Myers et al. (2007) 


1.47 


2 42+0-20 
^■^^-0.21 




[4] 


2QZ 


Porciani, Magliocchetti & Norberg (20041 




1.44 ±0.22 


0.63 ±0.10 




initial 






1.37 ±0.22 


0.64 ±0.10 




homogeneous 




1.47 


1.50 ±0.37 
1.35 ±0.23 
1.30 ±0.44 
1.54 ±0.51 


0.61 ± 0.15 
0.67 ±0.12 
0.70 ±0.24 
0.67 ±0.12 


decribed in text 


uniform 
good 
low reddening 
high reddening 


present paper 



fixed 7 on the scales > lQh~^ Mpc and we have found (3 
and 6 using Eqs. (O and ([8]). The results are given in Ta- 
ble 3 that contains the obtained values of infall and bias 
parameters and results of other authors with indication of 
the method, used for their estimation. One of them (denoted 
as [1]), containing simultaneous det ermination o f P and cos- 
mological model, was proposed by iHovle et al.l (l2002l ) and 
inclu des fitting of th e observed ^(a, tt) by theo retical model 



from iPeeblej l|l980l ) (see also iRatcliffe et"al] (I199S)) with 



as 

the relation for infall velocity from lHale-SuttonI Jl990l). The 
impro vement of this method proposed in Ida Angela et al.l 
l|2005h (denoted as [2]) allows to break the degeneracy be- 
tween JIm and p and is based in the consideration of the 
quasar cluste r ing ev olution. Another approach proposed in 
ICroom et al.l l|2005l ) (denoted as [3]) suggests the solution 
of the quadratic equation in bq with the shape of ^(r) gov- 
erned by the underlying dark matter distribut i on describe d 
by the analytic equ a tion f rom iHamilton et al.l (Il99ll . ll995h . 
In Ida Angela et al.l (|2008l 1 and IRoss et all (| 2009l) the last 
approach was applied with theor etical dark matter power 
spectrum from lSmith et all (|2003l ). One more method which 



was used in iPorciani. Magliocchetti fc Norber3 (|2004l ) (de- 
noted as [4]) includes using the relation between projected 
2pCFs of the quasars and that of the matter. As distinct 
from a number of authors (aimed, in particular, at estima- 
tions of the dark matter halo masses) we concentrate to the 
comparison with the linear theory; i.e. our results on bq are 
completely based on the interval r > 10 Mpc, where 
there is no problems with the linear theory. 



For test we also provided the same calculations for 
redshift-space 2pCF, P and b for homogeneous, uniform, 
good, low reddening and high reddening samples. The scat- 
tering of the resulting redshift-space 2pCF parameters and 
infall parameters with the corresponding values of so for 50 
homogeneous samples are presented in Figs. ll4l and ll5l Their 
mean values are presented in Tables 2 and 3 together with 
results for other samples. For all our samples errors of the 
parameters are jackknifes, except the homogeneous where 
errors are rms. All these results agrees well (within 1.3cr) 
with the results for initial sample. 
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Figure 15. Correlation lengtiies (for fixed slope value) and infall 
parameters for 50 homogeneous samples. 



5 RANDOM VELOCITIES DISPERSION 
5.1 Technique 

In this Section we discuss the 'Finger of God' effect as 
a tool to estimate the random component of the peculiar ve- 
locity V of quasars with respect to the CMB frame, i.e. the 
dispersion (v^). We neglect virialization and assume that the 
random velocities do not correla te with mass density varia- 
tions (contrary to iKaiseJ |l983) effect); such assumption is 
quite reasonable in case of quasars. If the quasars in a pair 
can be considered as independent, the line-of-sight pairwise 
velocity dispersion is (w^) = 2{v^), where one-dimensional 
pairwise velocity w is the difference of the peculiar veloc- 
ities (above the Hubble flow) vi, V2 of quasars in a pair: 
w = vi — V2. The effect of random velocities typically dom- 
inates at smaller separations, where the gravitational infall 
is less pronounced. Here, we neglect the effect of the gravi- 
tational infall. 

There are two major obstacles. Firstly, as it was men- 
tioned above, the SDSS data are incomplete in the region 
of linear comoving projected separations a < Lq = 2 Mpc. 
Nevertheless one can try to use an information from this re- 
gion by testing different distributions. Here we assume that 
the selection effects due to the fibre collisions do not distort 
considerably the distribution of quasar pairs from our sam- 
ple with projected comoving separations a < Lo = 2 Mpc 
and the redshift 0.8 < z < 2.2. Another obstacle is the 
poor knowledge of the 2pCF form in this region. This is the 
main point of this Section: to study the model dependence 
of the results for (w^)^''^. In order to obtain (v'^)^^^, we need 
2pCF for all distances. Therefore, we used the power-law 
form of the real-space 2pCF within the comoving distance 
range 2 < r < 50 Mpc with the parameters defined in Sec. [3] 
and assumed a rough estimate for r ^ Lq: 



(ro/r-r, r > Lo, 

J7(ro/io)^(Lo/r)^i, r Lo- 



(10) 



According to Sec. |3] we used 7 — 1.85, ro = 5.8 Mpc 
and 71, rj will be varied (71 < 3 to avoid divergence for 
r — > 0) . In case of ?? = 1 2 pCF is represen t ed by a continuous 
function. On account of iHennawi et abl l|2006l ) results, the 
rough choice 71 = 0, = 10 may be considered as plausible. 
However, even rather wild (rj — 100) choices do not lead 



to a drastic differences of results. Namely, we saw that the 
results are rather stable (taking into account errors) under 
different choices of 71, rj. 

Formula (|10|) allows to estimate the number of pairs 
with separation r. Then we restricted ourselves with the 
region a < Lq and use relation 



V'(n) = 27r / daaS, 



H^ 



to derive distribution of pair numbers (minus the random 
background) in line-of-sight separations H. 

We considered two kinds of the pairwise velocity dis- 
tribution: Gaussian and exponential ( Ratcliffe et al.l Il998l : 
iHovle et all I2OO2I : iHawkins et aD l200a r The pairwise ve- 
locity w contributes to the change of the comoving line- 
of-sight distance H between the quasars in a pair as fol- 
lows: H \n + d\, where 5 = w/[HoK{z)] with K{z) = 
[r2M(l + z)^ + 1 - Qm]'^^^/{1 + z). This coefficient grows 
from 0.85 to 0.97 within the interval 0.8 < jz < 2.2, hence 
we approximate it as K{z) ~ A"(1.5) ~ 0.9. 

First we assumed the expo nential pairwise velocity dis- 
tribution iRatcliffe et al.l [19981 ). Then we have probability 
density ^{5) of random variable 5 as follows 



$(5) = -exp{-a\5\) , 



V2KH0 

{w2)l/2 ■ 



In case of Gaussian distribution of velocities we have 



V2^ 



: exp 



2cr2 J 



2\l/2 



KHo 



(11) 



(12) 



The convolution of with <J> yields the distribution of pair 
numbers as a function of H. We used this to fit the observed 
pair numbers in different bins with the projected separation 
a < Lo — 2 Mpc after randomized background subtraction. 
The randomized background is calculated as described in 
Sec. [22] 

5.2 Results and discussion 

Because the number of pairs is not large, there is a 
considerable variation of results for different binnings, so 
we do not consider the numerical estimates as very reli- 
able. However, we may argue that the model dependence 
on 71 and parameters is rather insignificant. For illus- 
tration, we present the results (Table 4) of fitting the pair 
numbers (after randomized background subtraction) on ac- 
count of (|lll) and (Table 5) on account of (|12p for two 
different binnings having minimal jackknife dispersion esti- 
mates for {v^)^^"^ = {w^ /2Y^'^ . Binning 1 involves numbers 
of pairs with line-of-sight separations Pi from five intervals 
[2,11], [11,20], [20,29], [29,38], [38,47] (Mpc). Binning 2 
deals with separations from three intervals [2,17], [17,32], 
[32, 47] (Mpc). The bar charts in Fig.[l6]and[l7]we show cor- 
responding numbers after the average background subtrac- 
tion. Here we prefer to avoid separations less than 2 Mpc. 
Also for illustration in Fig. 1161 and 1171 we present the results 
of fitting for 71 = 1.85 and different values of r]. As one can 
see from the results, the values of (ii^)^''^ do not vary con- 
siderably under rather different (even unrealistic) choices of 
unknown parameters 71: the variations are of the order 
of the jackknife dispersion estimate. 
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Table 4. The random velocity dispersion {v'^)^^^ (km/s) for 
exponential distribution. Upper values correspond to Binning 
1, lower - to Binning 2. 





71 = 


71 = 1.85 


71 = 2.9 


v = i 


620±140 


660±150 


740±160 




460 ± 180 


500±180 


580±190 


n = 10 


720±160 


740±170 


750±170 




550±190 


580±190 


590±190 


jj = 100 


750±165 


750±165 


750±165 




590±190 


590±190 


590±190 



Table 5. The random velocity dispersion (ti^)^''^ (km/s) for 
Gaussian distribution. Upper values correspond to Binning 1, 



to Binning 


2. 








71 = 


71 = 1.85 


71 = 2.9 


■n = i 


580±90 


600±90 


670±90 




490 ± 120 


535±125 


600±110 


77 = 10 


650±100 


670±100 


690±80 




560±125 


590±130 


600±120 


T) = 100 


670±90 


670±90 


690±95 




580±130 


eOOitlOO 


600±120 



] 7, = 1.85, 77=1, <v'>""=660 km/s 
>7,= 1.85, rj=W, <v'>'"=740 km/s _ 
j7,= l,85, ?7=100, <v'>'"=750 km/sl 



20 30 

n, Mpc 



□ 7, = 1.85, 77=1, <v^>'^'=500 km/s _ 
o 7, = 1.85, 77=10, <v'>'"=670 km/s 
a 7, = 1.85, rj=100, <v'>'"=670 km/^ 



-2 - 




20 30 

n, Mpc 



Figure 16. Line-of-sight distance distribution of pair numbers in 
case of Binning 2 with fitted distributions for different values of 
71, rj, for exponential (top) and Gaussian (bottom) dis- 

tributions. The rectangles show the background subtracted pair 
numbers, the symbols with errors show the results of fitting. 



In fact the values {v^)^^"^ in Tables 5, 6 do not represent 
real quasar velocity dispersion, but a superposition of two 
components: (w^) = «err + (w^)) where (v^)^''^ is the true 
velocity dispersion and due to the redshift measure- 

ment errors. The last value correspond s to intrinsic emissio n 
lines shifts in quasars pointed out by I Groom et al.l (|2005l ). 
as it is an estimation from emission-line or cross-correlation 
technique of redshift measurement caused by the impos- 
sibility of precise determination of the line centre due to 
its large width . For detailed investigatio n of this effect see 
[Richards et"al] (|2002l ): IShen et al.l l|2007^ . The redshift mea- 



□ 7, = 1.85, 77=1, <v=>'" = 500 km/s 
o7,= 1.85, ?)=10, <v">"'=580 km/s 
a7,= 1.85, ?7=100, <v">'"=590 km/s 



20 30 

n, Mpc 



E r, _ 



□ 7, = 1.85, 77=1, <v'>'^'=535 km/s _ 
«7,= 1.85, 77=10, <v'>"" = 590 km/s 
a7,= 1.85, 77=100, <v'>"'=600 km/s 



20 30 
n, Mpc 



Figure 17. The same as in Fig. 1161 in case of Binning 2. 



surement errors from the SDSS data base for our sample is 
Az = 0.0019, which corresponds to Verr ~ cAz/(l + z) = 
230 km/s. In this case we would obtain too large value of 

Once work on th e first version of the text had been 
finished, the paper bv lHewett fc WildI l|2010l ) was appeared 
with the error estimate Verr = 600 km/s for SDSS data. 
In this case z-errors dominate in {v If we take the re- 

sults from Table 5, which yield smaller jackknife dispersion 
estimate for {v'^)^^^, than the residual upper estimate of 
{v^y^^ < (690^ - 600^)^/^ = 340 km/s seems to be more 
adequate. 

The values of the pairwise velocity dispersion (w^) = 
2{w^) corresponding to Tables 4, 5 are comparable with 
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that mentioned in previous studies of quasar clustering. 
lOutram et al.l l|200ll ) used the value 400 km/s obtain e d from 
the Hubble Volume Simulations; Ida Angela et al.l (|2005l ) 
fixed the value 800 km/s and found that this value gives an 
adequate result for s < 10 Mpc region and noted that 
it is dominated by the rms pairwise redshift error 600 km/s ; 
the same value 800 km/s was used bv lda Angela et al. (|2008|) 
and t he close value 690 km/s was used by Croom et alj 
l|2005i ). The last value was chosen as a mean one of the 
range 630-750 km/s; this is a quadrature superposition of 
(i) the real pairwise velocity dispersion recalculated with 
redshift estimated from the galaxies pa irwise velocity dis- 
persion 500 km/s iHawkins et all bOOSl ). (ii) redshift mea- 
surement error obtained from repeat observations, and (iii) 
velocity error due to intrins ic emission lines shifts in QSOs 
l|Richards et al.ll2002h . Onlv lMountrichas et all (|2009l ') treat- 
ing the velocity dispersion as a free parameter estimated 
its value for a sample of quasars and luminous red galax- 
ies (LRG) as 620 km/s and 727 km/s (two values are the 
result of two different approaches) and noted that they are 
consistent with the value 728 km/s which is a quadrature su- 
perposition of the values 800 km/s for quasars and 300 km/s 
for LRG from previous studies. 



on pairs with the projected separations a < 2 Mpc, which 
are not complete and the number of pairs involved in treat- 
ment of the velocity dispersion is not too large. In view 
of these circumstances we consider our results on (w^)^^^ 
rather as motivation to further studies to observe larger 
sample of ~ 1 Mpc quasar pairs with accurate redshift mea- 
surements. 
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6 CONCLUSIONS 

We have performed a detailed analysis of 2pCF using 
the quasar sample from the SDSS DR7 data for objects with 
0.8 < z < 2.2. Our main conclusions are as follows. 

(i) Using our approach for estimation of the real-space 
2pCF ^(r) of quasars we confirm that it is fitted well with the 
power-low model in the distance range 1 < r < 35 Mpc 
with the correlation length of ro = 5.85 ±0.33 Mpc and 
the power-low slope 7 = 1.87 ± 0.07. The results are in a 
good agreement with the results of the other authors on the 
same distance interval. However, this comparison suggests 
that the power-low model seems to be not valid on larger 
intervals. 

(ii) The redshift-space 2pCF of quasars has a break and it 
may be approximated well with the double power-low model 
with So = 6.43 ± 0.63 h'^ Mpc and 7 = 1.21 ± 0.24 for 
1 < s < 10 Mpc and so = 7.37 ± 0.81 ft"^ Mpc and 

7 = 1.90 ± 0.24 on scales 10 < s < 35 h^'^ Mpc wfth the 
mean values 5.83 ± 0.47 h'^ Mpc and 7 = 1.49 ± 0.12 for 
both distance ranges together. 

(iii) For distances s > 10 Mpc the parameter de- 
scribing the large-scale infall to density inhomogeneities is 
(3 = 0.63 ± 0.10 wfth the linear bias to be 6 = 1.44 ± 0.22 
that marginally (within 2cr) agrees with the linear theory of 
cosmological perturbations. 

(iv) We point out that the value of the quasars velocity 
dispersion varies almost within errors for a wide range of 
model parameters that determine 2pCF for r < 2 Mpc. 

We also performed investigations of our sample and 
showed that we can neglect boundary effects and all inhomo- 
geneities of the sample caused by non-uniformity of the sky 
coverage by the spectroscopic survey, reddening and condi- 
tions of observations. 

It is worth to note, that the values of {uP'Y^^ from Ta- 
bles 4, 5 have been calculated on account of SDSS DR7 data 
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APPENDIX A: DISTANCES DETERMINATION 

We consider the spatially flat Universe with Friedmann- 
Robertson- Walker (FRW) metric 

ds^ = C^dt'^ - a{t) [dx^ + X (dS^ + cos^ Sda'^)] , (Al) 

where a and S are right ascention and declination respec- 
tively, a{to) = 1 at the present epoch to. In case of the fiat 
AGDM cosmological model {Qm + Qa ~ 1) for an object 
with redshift z we have 

y(z) = ^ r ^ t(z) 

Ho Jo ^QA + nuii + cr 

The distances in the redshift-space are defined as if 
the objects were at rest with respect to the GMB. 
In order to define the distance r^jv from a quasar 
with coordinates {tq,Xq,<^q,Sq} to a neighbouring object 
with coordinates itN,XN,aN ,SN} we proceed as follows 
(llvashchenko fc Zhdanov! [2O10h . 

(a) Gonsider the Fermi frame (Synge 1960; Misner et 
al. 1973) of the quasar q with the fiducial trajectory a;^(t) = 
{t,Xq,o:q,Sq} that may be interpreted as the trajectory of 
an observer, which is at rest with respect to the GMB. 

(b) Introduce the hypersurface of constant time in 
corresponding normal Fermi coordinates with the origin 
{tq,Xq, ctq, Sq} . In our casc this is the space-like hypersur- 
face t = tq, which is orthogonal to the fiducial trajectory at 
the origin. 

(c) Find the invariant (proper) distance RqN from the 
origin to the intersection of the neighbour world line with the 
hypersurface; this is uniquely determined along the space- 
like geodesic connecting two points. In a general case we 
must suppose that the objects are not too far from each 
other to avoid the conjugate points. 

In our case RqN is easily determined along the geodesic 
on the hypersurface t = tq between the points {tq,Xq, CfqjSq} 
and {tq,XN , ctN , Sn}', then the comoving distance rgjv is: 

TqN = (1 + Zq)RqN = \J Xq + - '^XqXN COs(l/;,]v), (A2) 

where tpqN is the angle between objects q, N on the sky. The 
comoving linear projected distance from g to TV is 

a = XN sm{tpqN), (A3) 
and the comoving line-of-sight distance from g to A'^ is 
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n = Xiv cos(^qjv) - Xq, (A4) 

The above formulae determine distances in the redshift- 
space. Transition to the real-space involves the velocities of 
quasars with respect to cosmological background. However, 
in case of the projected distances the effect of velocities is 

negligible and these distances obtained from observations 
may be used to derive parameters of the real-space 2pCF. 

This paper has been typeset from a T^JC/ MJ^X file prepared 
by the author. 



